Atomic Gases at Negative Kinetic Temperature 
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We show that thermahzation of the motion of atoms at negative temperature is possible in an 
optical lattice, for conditions that are feasible in current experiments. We present a method for 
reversibly inverting the temperature of a trapped gas. Moreover, a negative-temperature ensem- 
ble can be cooled (reducing \T\) by evaporation of the lowest-energy particles. This enables the 
attainment of the Bose- Einstein condensation phase transition at negative temperature. 
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Amplification of a macroscopic degree of freedom, such 
as a laser or maser field, or the motion of an object, is 
a phenomenon of enormous technological and scientific 
importance. In classical thermodynamics, amplification 
is not possible in thermal equilibrium at positive temper- 
ature. In contrast, due to the fluctuation-dissipation the- 
orems 0, 01 , amplification is a natural phenomenon, and 
fully consistent with the second law of thermodynamics, 
in a negative-temperature heat bath Negative 
temperatures were introduced by Purcell and Pound 
to describe spin systems, and further generalized theoret- 
ically by Ramsey and Klein 0| . Useful introductions 
to thermodynamics at negative temperature are given by 
Landau and Lifshitz and by Kittel and Kroemer 
Negative temperatures are used to describe distributions 
where the occupation probability of a quantum state in- 
creases with the energy of the state. The prerequisite for 
a negative-temperature ensemble to be in internal equi- 
librium is that it occupies a part of phase space where the 
entropy decreases with internal energy. An ensemble in 
this part of phase space cannot be in thermal equilibrium 
with any positive-temperature system. An increasing oc- 
cupation number (and a decreasing entropy) can only 
exist up to a certain maximum energy. Unlike the mag- 
netic energy of spins, kinetic energy is not bounded from 
above in free space, which apparently rules out a nega- 
tive kinetic temperature. Here, we show that the band 
structure of a moderately deep optical lattice defines an 
upper bound to the kinetic energy of atoms in the first 
band. Therefore, atoms in an optical lattice form a sim- 
ple system in which thermahzation of the kinetic energy 
is possible at negative temperatures. 

Optical lattices, which are periodic potentials created 
through optical standing waves |l2l| , impose a band struc- 
ture on the motion of ultracold atoms. The kinetic energy 
of atoms in the lattice depends on their effective mass or 
band-mass, which, in analogy to the case of electrons 
in a semiconductor, can be strongly different from the 
normal rest mass. The effective mass can even become 
negative, as has been demonstrated in one-dimensional 
optical lattices 0, Q| . An isotropic negative effective 
mass can be realized in the vicinity of a band gap in 



a three-dimensional (3D) optical lattice. The bandgap 
is an energy range inside which no single-particle eigen- 
states exist. If the atoms are constrained to energies just 
below the bandgap the effective mass will be negative in 
all directions and the entropy of the ensemble will de- 
crease with energy. This is the condition for existence of 
a thermal ensemble with negative temperature j3| ■ 

For thermahzation to take place, the characteristic 
timescale for exchange of energy between particles must 
be shorter than timescales associated to heating and loss. 
In other terms, we need a high rate of intraband scatter- 
ing, with negligible interband scattering. Interband scat- 
tering is a collision process between two atoms in the first 
( "valence" ) band which causes one atom to be promoted 
to the second ("conduction") band while the "valence" 
band atom loses energy. Additional losses may be caused 
by interband tunneling, a one-body process in which an 
atom tunnels through a classically forbidden region, to 
re-appear in a different band. 

In this Letter, we show that all the requirements for 
thermahzation at negative temperature can be met in a 
3D optical lattice. First, we show that an optical lattice 
can have a 3D bandgap for atoms, at experimentally fea- 
sible lattice intensities. Moreover, interband scattering 
can be completely suppressed, and confinement of the 
negative-mass atoms can be easily accomplished using 
a magnetic trap for high-field seeking states. We calcu- 
late rates of interband tunneling and find they can be 
made negligible. Finally, we show how a trapped gas 
at positive temperature can be reversibly converted to a 
negative temperature. 

The motion of atoms in a simple cubic 3D optical lat- 
tice with a harmonic confining potential is described by 
the Schrodinger equation (See e.g., [l^ll5| ). 
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with Uut — t/o(cos^ kx + cos^ ky + cos^ kz). 

Here m is the bare atomic mass, a is the strength of 
the harmonic potential, Uq is the depth of the standing 
wave potentials that form the lattice, and k ~ 27r/A, 
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FIG. 1: Band structure of the one-dimensional Mathieu equa- 
tion. White area: Bands, grey area: gaps. Dotted hue: 
q = 0.559, dashed hne: q — 0.936, see text. 
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FIG. 2; "Valence" and "conduction" bands of a ID optical 
lattice a,t q — 0.936. Grey area; Gaps, white area; Bands. 



where A is the wavelength of the lattice light. The 
above Schrodinger equation is trivially separable by set- 
ting E = Ex + Ey + Ez- We now transform the equation 
to dimensionless units, 

x,y,z = kx,ky,kz (2) 
a = a/{Erk^) (3) 

a. = {Ex-^Uo)/Er (4) 
q = Uo/AEr, (5) 

with the lattice recoil energy Er = h?k'^ /2m. The one- 
dimensional equation of motion takes the shape of a mod- 
ified Mathieu equation, 

-— + {a^ - 2qcos2x)-^ + ax'^'i/ = 0. (6) 
dx'^ 

In the absence of a parabolic potential {a — 0), Eq. (|SJ) is 
exactly the Mathieu equation. Its solutions are the Math- 
ieu functions, which are propagating Bloch-Flouquet 
waves for values of a and q in the bands, and evanes- 
cent functions in the bandgaps. The relevant bands are 
depicted in Fig. ^ Adopting the usual convention 0| 
we designate the eigenvalues corresponding to the band 
bottoms by ao(g), ai{q), etc., and the band tops by bi{q), 
etc. (see Fig. If < \a\ ^ q/x, the lattice is locally 
only weakly perturbed by the parabolic potential, and 
the band structure will remain essentially intact. The 
energy positions of the bands become position depen- 
dent, in the same way as in a semiconductor in an ap- 
plied electric field. A negative value of a corresponds to 
a confining potential for negative-effective mass atoms, 
see Fig. |2 The area under the "valence" band in Fig. [21 
is classically forbidden as the kinetic energy of the atoms 
in this region is less than zero. The area between the 
bands is Bragg forbidden, i.e., no propagating modes can 
exist here due to interference. 



For the negative-temperature phase to be metastable, 
there must be a bandgap between the top of the "va- 
lence" band, where the negative mass states exist, and 
the bottom of the "conduction" band, where the effec- 
tive mass is positive. In a ID system (with the motion 
in the other two directions confined) this bandgap ex- 
ists for |g| > 0. However, in a 3D optical lattice, the 
bands may overlap. The maximum energy an atom in 
the "valence" band of a 3D lattice can have is 3bi{q), 
whereas the minimum energy in the "conduction" band 
is 2ao{q) + ai(q). The criterium for the existence of a 
bandgap is 36i(g) < 2ao{q) + ai{q), which is fulfilled for 
q > 0.559. 

If the 3D bandgap is narrow, collisions between two 
"valence" band atoms can promote an atom to the "con- 
duction" band, where it has positive effective mass. This 
interband scattering will not only lead to loss of atoms, 
but also to loss of energy, which is equivalent to heat- 
ing (increasing |T|). We show now that interband scat- 
tering becomes energetically forbidden in two-body pro- 
cesses at high enough q. The maximum energy of the 
initial state, consisting of two "valence" band atoms, is 
6bi{q). The minimum energy of the final state, consist- 
ing of one "valence" and one "conduction" band atom, is 
[3ao(g)] -|- [2ao(g) +ai{q)]. The interband scattering pro- 
cess is energetically forbidden if 65i(g) < 5ao(g) -I- ai(g), 
which is the case iiq > 0.936. This value of q corresponds 
to a lattice depth almost an order of magnitude less than 
the depth employed to create a Mott insulator phase |0| . 
In such a moderately deep lattice, thermalization rates 
and inelastic collision rates are all of the same order of 
magnitude as the respective rates at positive temperature 
in the absence of a lattice 0| . Through elastic collisions, 
the ensemble will always relax to a distribution with the 
highest possible entropy at its given energy, in this case 
a negative-temperature distribution. 

To investigate how the gas thermalizes in the negative- 





FIG. 3: Density of states p{E) in a three-dimensional lattice 
with a = —10"'^. The energy E is taken with respect to the 
top of the first ("valence") band. 



temperature phase space, it is useful to consider the semi- 
classical density of states (DOS). The DOS can usually 
be found by taking the energy derivative of the Wentzel- 
Kramers-Brillouin (WKB) phase integral, a procedure re- 
cently applied to a one-dimensional optical lattice with 
harmonic confinement (positive a) |l8|. For negative a 
the WKB method cannot be used straightforwardly to 
calculate wavef unctions, as the "Bragg reflection" bound- 
ary condition at the outer turning points behaves differ- 
ently from classical turning points. The WKB formula 
for the DOS does not depend on the boundary conditions 
and can be straightforwardly calculated in the one-band 
model of .18]- Within the one-band and WKB approxi- 
mations, the DOS Pa at negative a is related to the DOS 
calculated for positive a by 



(7) 



where in the case of negative a, the top of the valence 
band is taken as the energy zero. The 3D DOS crosses 
over smoothly from a quadratic behavior at small neg- 
ative energy to approximately square-root behavior at 
larger negative energy, as shown in Fig. |21 The DOS is 
nonzero for any energy lower than the band edge at the 
trap center, so in principle the negative temperatures are 
not bounded in an infinite lattice. In fact, this asymmet- 
ric behavior of the DOS makes thermal equilibrium at 
a positive temperature impossible, which is in contrast 
to spin systems 0, where thermal equilibrium at posi- 
tive and negative temperature can be studied in the same 
field configuration. 

We now examine the escape of the atoms through in- 
terband (Zener) tunneling. The trapped negative-mass 
"valence" band states are mixed with anti-trapped "con- 
duction" band states at the same energy, as the evanes- 
cent tails of these states extend through the Bragg for- 
bidden region. The Zener tunneling rate Fzcnor may be 



FIG. 4; Method for reversibly transforming a positive- 
temperature gas (a) into a negative-temperature gas (f). Each 
step is reversible, see text, (a) trapped gas at positive T. 
(b) a strong optical lattice, (c) a strong optical lattice with- 
out trap, (d) a strong optical lattice with changed asc, (e) a 
strong optical lattice with negative trap, (f) an optical lattice 
at g = 0.936 with a negative-temperature gas. 
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where kB\och{x) is the Bloch wavevector of the Mathieu 
equation corresponding to the energy E — ax^ , and wa is 
an attempt frequency of the order of the trap oscillation 
frequency (typically < 1 kHz). The resulting tunnel- 
ing time constant is of order 10^ seconds at a = 10^"^. 
For 10"'' < a < 10~^ we checked Eq. |H1 numerically by 
high-precision quadrature of Eq. El and found agreement 
within an order of magnitude. 

The conditions for existence of a negative-mass Bose- 
Einstein condensate (BEG) were recently studied by Pu 
and coworkers 0|, they also investigated the possibility 
of creating such a condensate by phase imprinting, and 
found this likely to be an inefficient mechanism as many 
atoms are transferred to other bands. Other mechanisms, 
such as accelerating the atoms uniformly to complete 
one-half Bloch oscillation, may be more efficient, espe- 
cially since in a BEG nonlinearity due to the interactions 
counteracts dephasing |^. Using these methods, 
the negative-mass BEG can be created in a mechanically 
unstable state as the interatomic scattering length age 
remains positive [Tsf . 

We propose a radically different mechanism for the pro- 
duction of negative-temperature gases, based on the re- 
versible creation of a Mott-insulator phase from a BEG 
|17|. The idea is to "freeze" a positive temperature gas 
into a Mott insulator, by turning on a strong optical lat- 
tice. The Mott phase is then "molten" in a different trap 
configuration and with different interatomic interactions, 
to produce a negative-temperature gas. The procedure 
is outlined in Fig. ^ we will describe it here for a pure 
BEG but there is no a priori reason why it would not 
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work when some thermal excitations are present. We 
start with a trapped BEC at a > and a^c > (a), and 
turn on a strong optical lattice, as is done in ^ Mott 
insulator phase will form in which the atoms are confined 
to their local potential wells (b). We assume one atom 
is present per lattice site. Now we remove the trap (c), 
which can be done rapidly without creation of entropy, 
as the atoms are tightly bound to their respective lattice 
sites. Then we change the sign of a^c (c,d), for example, 
by tuning the magnetic bias field near a Feshbach reso- 
nance (See e.g., |23). Since the atoms are at different 
lattice sites and do not interact, this can be done rapidly 
without causing losses. The system passes through the 
point where both a^c and a are zero, here the temper- 
ature of the Mott phase is not defined, but the entropy 
is. This passage transforms the Mott insulator from the 
lowest-energy many-particle state to the highest energy 
in the "valence" band. The confining potential is recon- 
figured to change the sign of a (e). This can be done 
either by changing the confinement field, or by chang- 
ing the internal state of the atoms. A magnetic field 
minimum will trap negative-mass atoms in their lowest 
Zeeman state [l^. No change in the spatial distribution 
of atoms will take place due to the extremely low tunnel- 
ing rate in the Mott insulator, so this step can proceed 
rapidly. When the lattice is slowly reduced in strength 
(f), the Mott insulator will melt into high energy states 
in the "valence" band, i.e., a negative-mass, negative- 
temperature gas. This approach seems practicable for all 
ultracold bosonic atoms which have suitably broad Fes- 
hbach resonances, such as ''Li, ^^Na, *^Rb, and ^"^^Cs. A 
detailed analysis of the reversibility criteria for each step 
will be published elsewhere. 

The atomic gas can be cooled by forced evaporation. 
In the case of a magnetic trap, the evaporation can pro- 
ceed through radio frequency (RF) induced spin flips to 
an untrapped (magnetic quantum number to = 0) state. 
One should tune the RF to selectively remove atoms with 
lower-than-average energy, so that the average energy per 
atom will increase. This leads to an increase of T, and a 
reduction of the entropy per particle S/N . A reduction 
of the entropy per particle is equivalent to an increase 
of the phase space density. Therefore we use the term 
cooling for any process that lowers |r|. At low enough 
|r| the system can be described by a (negative) power 
law density of states, and the theory of evaporative cool- 
ing in a power law density of states can be applied. 
(One should take care to reverse the signs of to, E and T 
in Ref.'23'l wherever applicable.) The condition for BEC 
in a 3D system is nA^ — 2.61, where n is the number 
density of atoms and A^ — 2TTh?m~^k^^T^^ , where ks 
is Boltzmann's constant. Since both the effective mass 
TOeff and T are negative, A is a real number. Achieving 
the BEC phase transition at negative temperature seems 
possible, either by heating (extracting energy), starting 
from a pure BEC, or by cooling, starting from a negative- 



thermal gas. 

The negative value of asc, which is necessary to create 
the negative-mass BEC from a Mott insulator, is also a 
condition of stability for the negative-mass BEC. In fact, 
the negative-mass BEC is likely to collapse if Osc is posi- 
tive 15]. It is very interesting to consider the possibility 
of superfluidity at negative a^c in a negative-temperature 
system. One expects the superfluid phase to allow for dis- 
sipationless motion of objects, unless a critical velocity 
is exceeded, in which case negative dissipation sets in. 
Further properties of negative-temperature superfluidity 
are hitherto unexplored theoretically as well as experi- 
mentally. 

In conclusion, in an optical lattice combined with a 
parabolic potential, a large phase space in which atoms 
can thermalize to negative kinetic temperatures is open 
for exploration. 

I have greatly benefited from inspiring discussions 
with Claude Cohen- Tannoudji, Ad Lagendijk and Willem 
Vos. 
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